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FINITE ELEMENT QUASI-INTERPOLATION AND BEST 
APPROXIMATION* 

ALEXANDRE ERNt AND JEAN-LUC GUERMOND* 


Abstract. This paper introduces a quasi-interpolation operator for scalar- and vector-valued 
finite element spaces constructed on affine, shape-regular meshes with some continuity across mesh 
interfaces. This operator gives optimal estimates of the best approximation error in any L^-norm 
assuming regularity in the fractional Sobolev spaces where p G [l,oo] and the smoothness 

index r can be arbitrarily close to zero. The operator is stable in leaves the corresponding finite 
element space point-wise invariant, and can be modified to handle homogeneous boundary conditions. 
The theory is illustrated on /f{curl)- and /f{div)-conforming spaces. 
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1. Introduction. Consider a shape-regular sequence of affine meshes {Th)h>o 
approximating a bounded polyhedral Lipschitz domain D in and a sequence 
of finite element spaces {P{Xh))h>o based on this mesh sequence, composed of ei¬ 
ther scalar- or vector-valued functions, and conforming in some functional space W 
where some continuity across mesh interfaces is enforced. Estimates of the best- 
approximation error in P{Th) are invariably invoked in the convergence analysis of 
finite element approximations. When the exact solution is smooth enough, the canon¬ 
ical interpolation operator in P{7h) can be used to obtain decay estimates of the 
best-approximation error in terms of the mesh-size. However, in many practical situ¬ 
ations, the exact solution only sits in a Sobolev space (for some p G [1, oo]) 

where the smoothness index r can be close to zero. In this case, an alternative quasi¬ 
interpolation operator must be invoked to estimate the best-approximation error. 

The aim of this paper is the construction of a quasi-interpolation operator that 
is stable in leaves P{Th) point-wise invariant, and approximates (quasi-)locally 
and optimally functions in W'^’P{D), for all p G [1, oo] and r arbitrarily close to zero. 
Moreover, the construction can be modified to enforce homogeneous boundary con¬ 
ditions that are legitimate in W. The main examples we have in mind are R^-, 
ff(curl)-, and i?(div)-conforming finite element spaces. Let us emphasize that, for 
vector-valued elements, the present construction is not a substitute to the notion of 
commuting bounded cochain projection introduced in the framework of Finite Ele¬ 
ment Exterior Calculus, but instead a vital complement to it. Indeed, as shown in 
Arnold et al. Theorem 5.6], the approximation error for the commuting projec¬ 
tion is bounded by the best-approximation error in P{7h), up to a uniform constant. 
Therefore, the present quasi-interpolation operator shows that commuting bounded 
cochain projections converge optimally for low-regularity solutions. Another impor¬ 
tant consequence of the present construction in the case of vector-valued elements is 
that the decay rates of the best-approximation error only involve the solution semi- 
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norm in without the need to invoke a bound on the curl or the divergence 

of the solution. 

Let us put our results in perspective with the literature. The construction of quasi¬ 
interpolation operators with the above properties for iL^-conforming finite elements 
is a well-studied field. In two space dimensions on triangular meshes, we mention the 
early work of Clement Q; see also Bernard! and Girault Q where the construction 
is modified so that the resulting operator leaves the finite element space point-wise 
invariant. In a landmark paper by Scott and Zhang (^ . an alternative construction is 
proposed in any space dimension for iJ^-conforming finite elements with and without 
homogeneous boundary conditions; see also Girault and Lions 15| (and [H, p. 491]) 
for a modification leading to L^-stability. The so-called Scott-Zhang interpolation 
operator leads to optimal decay estimates of the approximation error for functions in 
The extension of this result to functions in W^’P{D) with r G (i, 1) has been 
studied only recently by Ciarlet Q using the original Scott-Zhang operator. To the 
authors’ best knowledge, in contrast to iL^-conforming elements, quasi-interpolation 
results for iT(curl)- and iT(div)-conforming elements are missing in the literature. 

The main ingredient for the Clement and Scott-Zhang construction is a regulariza¬ 
tion of functions based on macroelements consisting of patches of elements. We adopt 
here a somewhat different route. The key idea is a two-step construction, namely a 
projection onto the (fully discontinuous) broken finite element space followed by an 
averaging operator that stitches the result in the spirit of Oswald [2^, Eqs. (25)-(26)]. 
This way, patches of mesh cells are only involved when handling discrete objects, 
thereby avoiding the delicate question of reference patches frequently used in the lit¬ 
erature. Averaging operators are often invoked in the literature in various contexts; 
we mention, in particular a posteriori error estimates [ll[l^[l3|. preconditioners 22). 
stabilization tec hniq ues [ 3 , and discontinuous Galerkin methods with improved sta¬ 
bility properties Ifl, 1^ • The novel feature here is the analysis of averaging operators 
in the vector-valued case and the application to quasi-interpolation. We emphasize 
that we devise a unifying framework for the analysis that works irrespective of the 
nature of the degrees of freedom, i.e., whether they are nodal values, integrals over 
edges, faces, or cells. Essentially all the finite elements that are traditionally used 
to build iT(curl)-, and JT(div)-conforming finite element spaces match the few 
assumptions of the unified analysis. 

The paper is organized as follows. In we introduce the notation and construct 
a sequence of abstract finite element spaces conforming in some functional space W. 
Key assumptions are identified and listed. A local interpolation operator stable in 
is constructed in ]J3) An abstract averaging operator acting only on discrete functions 
is introduced and analyzed in (21 The final quasi-interpolation operator is constructed 
in 21 without enforcing any boundary condition. The question of enforcing boundary 
conditions is addressed in 21 Finally, 21 contains two technical results on fractional 
Sobolev spaces that are of independent interest: a Poincare inequality and a trace 
inequality. 


2. Finite elements. In this section we introduce some notation and construct a 
sequence of abstract finite element spaces, conforming in some functional space W. In 
the entire paper the space dimension is denoted by d and the domain D is a bounded 
Lipschitz polyhedron in 

2.1. Meshes. Let {Th)h>o be a mesh sequence that we assume to be affine and 
shape-regular in the sense of Giarlet. We also assume for the sake of simplicity that 
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the meshes cover D exactly, that they are composed of convex cells, and that they 
are matching, i.e., for all cells K, K' G Th such that K ^ K' and K n K' ^ the set 
K n K' is a common vertex, edge, or face of both K and K' (with obvious extensions 
in higher space dimensions). By convention, given a mesh Th, the elements K G Th 
are closed sets in 

We assume that there is a reference element K such that for any mesh Th and any 
cell K G Th, there is an affine bijective map between K and K, which we henceforth 
denote Tr- : K ^ K. Since Tk is affine and bijective, there is an invertible matrix 
Jif G such that 

Tk{x)-T xiy) = Ik{x - y), Vx,y G K. (2.1) 

In what follows, we denote points in R*^, R^-valued functions and R'^-valued maps in 
boldface type, and we denote the Euclidean norm in R'^ by ||•||^ 2 (]g<^), or Il'Ufa when 
the context is unambiguous. We abuse the notation by using the same symbol for 
the induced matrix norm. Owing to the shape-regularity assumption of the mesh 
sequence, there are uniform constants c^, such that 

|det(JK)| = \\rK^h^<c^h-^\ (2.2) 

where Hk is the diameter of K. Recall that = — and = —h^ for meshes 

Pk PK K 

composed of simplices, where px is the diameter of the largest ball that can be 
inscribed in K, is the diameter of K, and is the diameter of the largest ball 
that can be inscribed in K. 

2.2. Finite element generation. We are going to consider various approxima¬ 
tion spaces based on the mesh sequence {Th)h>o- Again for the sake of simplicity, 
we assume that each approximation space is constructed from a fixed reference fi¬ 
nite element (K, P, E). The reference degrees of freedom are denoted {cti, ..., 
and the associated reference shape functions are denoted {0i,..., by definition 
di{0j) = 6ij, \/i,j G {lirish}- We denote N := {lirish} to alleviate the notation. The 
shape functions are R'J-valued for some integer q > 1. We henceforth assume that 
P C 1T^’°°(A'; R'J) (recall that P is typically a space of polynomial functions, but we 
do not require this assumption here). 

We assume that there exists a Banach space V{K) C such that the 

linear forms {cti, ... ,(?nsh} oan be extended to £(F(A');R), i.e., V{K) is the domain 
of the degrees of freedom; see [3, P- 39]. Then, we define : V{K) -G P, the 
interpolation operator associated with the reference finite element {K, P,E), by 

Vx gK, TvG V{K). (2.3) 

i&M 

By construction, G C{V{K)\ P), and P is point-wise invariant under 

We now address the question of constructing finite elements for any mesh cell 
K G Th- We assume that there exists a Banach space V{K) C L^{K\R'^) and a 
bounded, bijective, linear map between V{K) and V{K): 

il^K ■■ V{K) 3v ^ %I}k(.v) G V{k). (2.4) 


We then set 


Pk ■■= {p = i^kip) \ p&P}, 

Er- := {crK.i}iejG s.t. aK,i = o tpx- 


(2.5a) 

(2.5b) 
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Proposition 2.1 (Finite element). The triple (K, Pk ^Tjk) is a finite element. 

Proof. Note first that dim(PR') = dim(P) = righ since ipK is bijective. Moreover, a 
function p G Pk such that aK,i [p) = 0 for alH G TV is such that ifK [p) = 0 by the uni- 
solvence property of the reference finite element; hence, p = 0. Finally, the linear forms 
aK,^ are in £(P(iF);R) since \crK,i{v)\ < \\'^^\\civ{Ky,R)\\^K\\c(v{Ky,v{K))Mv{K), for 
all V G V{K). □ 

The above definitions lead us to consider the canonical interpolation operator 
associated with the finite element {K, Pk, T,k). 


Ik{v){x) = ^aK,i{v)9K,i{x), Vx € K, Vv€V{K), (2.6) 

ieAf 

where we have set 9k, i ■= '4’~f^{9i). Note that Ik G C{V{K)-,Pk) and that Pk is 
point-wise invariant under Ik- 

Since the mesh is affine, we assume that ipK has a simple structure; more precisely, 
we assume that there is a, qxq invertible matrix Ak such that 


'ipK{v) = Ak{v o Tk), 


(2.7) 


which implies that i^k can be extended as a map from L^{K]W>) to L^(iF;R'*). The 
following classical result shows that fiK maps M'^) to for all / G N 

and all p G [l,oo]. In particular, this implies that Pk C IT^’°“(iF; R.'*). We use the 

notation \V’K\mYi,p(^K-,Ri)-,w‘’P(K-,Ri)) ®'^P«ewcp(R';R<i) similar 

notation for \'f’K^\c(^\Y‘-p{K-Ri)-w‘’P{K-Ri))^ where it henceforth implicitly understood 


the denominator is not zero. 

Lemma 2.2 (Bound in Sobolev norms). Let I G N. There is a uniform constant 
c depending on the shape-regularity of the mesh sequence {Ph)h>o on I such that 
the following holds: 


\'^^\c{W‘'P{K;Ri);W‘'P{K;Ri)) — cllAKH^a ||Jr -||^2 |det(J_R-)| P, 

\'^K l£(WCP(X;R9);Wi,P(i<';R9)) — c||A^ ||f2 ||Jjf ||^2 |det(jR-)|p, 

for all K (yTh and all p G [1, oo] (with z^p = 1, Vz > 0 if p = oo). 

Proof. For any multilinear map A G Adj(R'^,... ,R'^;R'^), let us set 


II^IIaIi(R‘‘....,R‘^;R‘') ■ — 


sup 

iyi,....yi)&R‘‘x...xR‘^ 


l | A ( t / i ,..., I /;)||^2 

'\yi\\t^ ■ ■ ■ WVlWl^ ' 


(2.8a) 

(2.8b) 


Then, denoting by D^ifK{v) the /-order Frechet derivative of ipK at v, the assump¬ 
tion (UTl) implies that ||TlVR:(^^)ll 7 Wi(R‘i....,R‘i;R‘i) < l|Aif||f 2 ||i:)'(i;oTK)||^^(jjd_^^^_jjd.R,) 
for all / G N. Then, standard results about the transformation of seminorms in 
the Sobolev space W^’’^ using the pullback by Tk lead to (12.81) . see e.g., Ciarlet [ 1 , 
Thm. 3.1.2] or Ern and Guermond [l^. Lemma 1.101]. □ 

Corollary 2.3 (Bound on Ak). Assume that there is a uniform constant c so 

that 


(2.9) 
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Then the following holds for all s, m G N, all K G Th and all p G [1, oo]; 

\'^K \c{W"'’P{K)-W’^’P{K))\'^k\c{W^’P{K)-W^’P{K)) — ' ( 2 - 10 ) 


Proof Combine (|2.8I) with (12.21) and (|2.9D . □ 

Lemma 2.4 (Norm equivalence). There exists a uniform constant c such that 

lk||L=o(^.R,) < c||A^^||f2 |crA.*(w)|, (2.11) 

iGA /" 

Y, \’^K,i{v)\ < c||AK||£2||u||ioo(;f.R5), (2-12) 

iGA / 

for all V G Pk and all K €Th. 

Proof. We prove (12.1111 : the proof for the other bound is similar. We observe that 
||w||l“(A;R<!) = \\'^k^{'^k{v))\\lpp{K-Ri) < 

< c||A^1||^2 Y \M'^Kiv))\ = c||A^1||^2 Y \^K,i{v)V 

ieAf ieAf 

using the definition of ipK, norm equivalence in P, and (I2.5bl) . □ 

2.3. Abstract finite element spaces. Let {{K,PK,T,K)}KGTh be a 7/t-based 
family of finite elements constructed as in Proposition 12.11 We consider the broken 
Sobolev spaces W1 -p(7;;R«) := {v G Lp(P;R«) | vik G \/K g Th}, 

p € [1, oo], and introduce the broken finite element space 

P'^iTh) = {vh G Li(P;M«) I vuiK € Pk, VA G %]■ (2.13) 

Since Pk C ^^’“(A;R?), we infer that P'^{Th) C Wi’“(7;;R«). 

Our aim is to define a subspace of P'°{Th) by means of some zero-jump condition 
across the interfaces separating the elements. We say that a subset F C P is an 
interface if it has positive (d—l)-dimensional measure and if there are distinct mesh 
cells Ki,Kr G Th such that F = dKi (~l dKr- The numbering of the two mesh cells is 
arbitrary, but kept fixed once and for all, and we let np he the unit normal vector to F 
pointing from A/ to A^. This defines a global orientation of the interfaces. We denote 
by riKi and riKr the outward unit normal of A; and Kr, i.e., np = ripi = —npr- 
We say that a subset P C P is a boundary face if it has positive (d—l)-dimensional 
measure and if there is a mesh cell K € Th such that F = OK id dD, and we let np he 
the unit normal vector to F pointing outward P. The interfaces are collected in the 
set P)(, the boundary faces in the set pf, and we let = P)( U pf. We now define 
a notion of jump across the interfaces. Recalling that functions in W^'^{Th','^'^) have 
traces in L^(clA;R'J) for all A G Th, let F G Tf be a mesh interface, and let A;, A^ 
be the two cells such that F = dKi (~l dKp, the jump of G W^'^{Th', R"^) across F is 
defined to be 


| u 1 f (£ c ) = v\Ki{x) - v\K^{x) ee X e F. (2.14) 

In what follows, we need to consider the jump of only some components of v across 
P; we formalize this by introducing bounded linear operators '^k,f ■ W^’^(A;R'J) —>■ 
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L^(F;]R‘), for all K G Th: all face F £ that is a subset of dK, and some t > 1, as 
follows: 


WF(a=) = lKi,Fiv\Kt){x) - -fKr,Fiv\K^){x) Si X £ F. (2.15) 

We assume that ||f]J'(3:)| < ||f]F(a5)|, sea: G F, for all v £ Since 

functions in have zero jump across interfaces (see, e.g., [ll|, Lemma 1.23]), 

we infer that 


V £ = 0, VF G F^ (2.16) 

With this setting, we introduce the 

P{Th) = {vh £ P^{%) I = 0, VF G Fn- (2.17) 

2.4. Finite element examples. The present theory is quite general and covers 
a large class of scalar- and vector-valued finite elements. For instance, it covers finite 
elements of Lagrange, Nedelec, and Raviart-Thomas type. To remain general, we 
denote the three reference elements corresponding to these three classes as follows: 
(£, pg, pg), {K, P\E^) and (£, E^). We think of {K, Ps, ps) as a scalar-valued 

finite element (g = 1) that has some degrees of freedom which require point evaluation, 
for instance (FT, F®, F®) could be a Lagrange element. We assume that the finite 
element {K,P^,E^) is vector-valued {q = d) and some of its degrees of freedom 
require to evaluate integrals over edges. Typically, {K,P^,E^) is a Nedelec-type 
or edge element. Likewise, the finite element {K, P^, E^) is assumed to be vector¬ 
valued (g = d) and some of its degrees of freedom are assumed to require evaluation 
of integrals over faces. Typically, {K,P‘^,E‘^) is a Raviart-Thomas-type element. 
It is not necessary to know the exact nature of the element that we are handling 
at the moment. We denote by V^{K), V^{K), V^{K) admissible domains of the 
degrees of freedom in the three cases. Let p £ [l,oo]. The above assumptions imply 
that we can choose V^iK) = W^’P{K) with s > V^K) = W^’P{K) with s > 
and V'^{K) = W‘^’P{K) with s > F Actually, when p = 1 we can choose 
V^{K) = W‘^’^(Fr) (since ^ C^{K), see e.g., Ponce and Van Schaftingen 

(iljjl. V‘^{K) = (since functions in have a trace in F^(9F')), and 

V^{K) = (since functions in W^“^’^(Fr) have traces in on the one¬ 

dimensional edges of J’f). 

Let 7/t be a mesh in the sequence {Th)h>o and let FT be a cell in Th- We denote 
by Pi-, Ppi Pa the linear maps that are classically used to generate the above 
finite elements, i.e., is the pullback by Tk, and ipK ipR 3-'^® the contravariant 
and covariant Piola transformations, respectively. All of these maps fit the general 
form (12.7L i.e., 

A^ = l, i^%iv)=voTK, (2.18a) 

A5f=jP, ii'k{v)=ll{voTK), (2.18b) 

Ap = det(JF') = deti^K) oTk). (2.18c) 

Note that c = 1 in (12.91) for the above examples. 
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The corresponding broken finite element spaces are: 

pg.b( 7 -^) ^ g I g ps, viC e r4, (2.19a) 

p-’^{Th) = {Vh e l\d) I ^‘k{^h\K) e P\ Vif e r4, (2.i9b) 

P^'^iTh) = {vh e I ^%{vh\K) e V/^ G Th}. (2.19c) 

This leads us to consider the following 7 -traces: 

7l:,F(^|if)(®) = '^\k{x), Va: G F, (2.20a) 

1 k.f{'^\k)[x) = v\K{x)xnF, Va; G F, (2.20b) 

1k.f{'^\k){x) = v\K{x)-nF: Va: G F, (2.20c) 

and the following conforming finite element spaces: 

P^{Th) = {vh e F8’4r4 I Kll^ = 0 , VF G Fa, (2.21a) 

p"(rft) = {vh G p^’^ra I KFf = 0 , vf g Fa, ( 2 . 21 b) 

p‘'(ra = {vn G p'^’^ra I kIf = 0 , vf g Fa, ( 2 . 21 c) 


where we slightly simplified the notation by using instead of etc. Upon 

introducing the spaces U® := {u G L^{D) \ Vv G L^{D)}, := {i; G L^{D) | Vxu G 

L^{D)}, V‘^ := {^; G L^D) \ V-v G L^D)}, we have PS(F) := P^’^^iTh) n U®, 
P^{rh) := P^’^^iTh) n Vp P‘^(r^) := P‘^’'"(r^) n that is to say, the finite element 
spaces P^{Th), P^{Th), P'^{Th) are conforming in U®, V''^, V''^, respectively. 

Let us introduce the canonical interpolation operators If^, such that 

^hi'^)\K = P%{vik), Phi'^)\K = ^P('W|f), ^hi'^)\K = 2:p(v|if). The considerations 
in fl2.4l show that it is legitimate to take 1U®’^(F), s > |, for the domain of X®, 
W^’P{D), s> for the domain of X^, and W®’P(F), s > i, for the domain of X^ 
i.e., the canonical interpolation operators X|, X)( and X)( are not stable in any LP{D) 
space (or Lp{D)). The objective of this paper is to construct quasi-interpolation op¬ 
erators mapping onto the spaces P^{Th), P^{Th) and P'^{Th) that are stable in L^{D) 
(or L^{D)) and have optimal approximation properties with and without boundary 
conditions. 

2.5. Summary of the assumptions. Let us now summarize the assumptions 
that will be used in the rest of the paper. Henceforth {Th)h>o is a shape-regular 
sequence of affine, matching meshes so that m and ( 1 ^ hold. We also assume 
that the map V'ic satisfies ( 1 ^ and (12.91) . {{K,PK,'PiK)}KGTh is a F-based sequence 
of finite elements constructed as in Proposition 12.11 In view of approximation, we 
define k to be the largest natural number such that [Pfc.d]"^ C F, where Fk,d is the real 
vector space of d-variate polynomials functions of degree at most k, and we assume 
that P C IU'=+i’°°(F;K9). 

We assume that we have at hand a notion of 7 -jump across mesh interfaces as 
described in (12.31 The finite element space P{Th) is the subspace of the broken finite 
element space P^{Th) characterized by zero 7 -jumps across interfaces, see (12.171) . 
Finally, two important assumptions relating the degrees of freedom to the 7 -jump 
and 7 -trace are the estimates (j4.2l) and (16.61) below. 

In what follows, c denotes a generic positive constant whose value may depend on 
the shape-regularity of the mesh sequence {7h)h>o and on the reference finite element 
{K, P, S). The value of this constant may vary from one occurrence to the other. 
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3. L^-stable local interpolation. In this section we extend the degrees of 
freedom in order to be able to approximate functions that are only integrable. 

3.1. Extension of degrees of freedom. Let us consider pi G P, i € Af, be 
such that 

-i- / {Pi,p)e'^{Ri)dx = ai{p), VpeP. (3.1) 

\K\ Jk 

Note that pi is well defined since it is the Riesz representative of ai in P when P is 
equipped with the L^-scalar product weighted by l/|if|. This leads us to define 

^t(p) ■= 7^7 / (Paw)^ 2(R9) dx, Vv G L^{K-,W). (3.2) 

\K\ Jk 

Note that the assumption P C L°°{K]W>) implies that pi G which in 

turn implies that all the extended degrees of freedom are indeed bounded 

over L\K-,R‘i) since pf < \K\-^\\pi\\^^^f^,^^y In passing we have also 

proved that < \K\-^\\MLp’iK-,Ri) forallp G [l,oo], where d + T = 1. 

We then define 

■= E S L\K-R'i). (3.3) 

i&M 

P is point-wise invariant under since a\ (p) = Oi (p) for all p G P and all i G A/”. 

Let K gTh and let {K, Pr-, S/^) be a finite element constructed as in (12.51) . Note 
that the assumption (12.711 implies that '0_R-(P^(Pr;]R‘?)) = P^(iL;R'^). We then extend 
the degrees of freedom in Tjk to L^(Pr;R'*) by setting 

(3.4) 

The above definition leads us to define 

:= Vn G pi(iL;R«). (3.5) 

i&M 

Proposition 3.1 (Stability, commutation, invariance), (i) There exists a uni¬ 
form constant c such that \\iTl^\\c{LP(K;m);LP(K-,Ri)) < c, for all p G [l,oo] and all 
K G Th', (ii) Tji commutes with ipK; (Hi) Pk is point-wise invariant under 1^^.. 

Proof. Using the triangle inequality in (13.51) . the fact that 6k, i ■= and 

i ■— (see (I3.4p i. and finally using (j2.8ll with I = 0 and the assumption (j2.9ll . 

we infer that 


l | T ^ IU ( LP ( A '; R ‘ i ); LP (_ fs '; R < i )) 

C{LP(K-,Rp)-LP{K-,R'i)) \\i’K^iOi)\\LP{K-,Rp) 


ieJ\f 

< It/'A 


C(,LP{K-,RP)-,LPiK: 


IV'. 


K \c{LP{K-,Ri)-,LP{K-,RP)) 


El 


b \IC{LP{K-,B 






< c ||Ai^||^2||A^ 11^2 ^ Iki I|£(LP(A;R 9 );R)II^*IIlp(A;R<i) 


iSTV 


< c||JJK||f2||J^^||f2|.P| ^ ^ IIPi|lLP'(A;R<l)- 

iGAf 
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The conclusion readily follows from the shape-regularity assumptions. To prove the 
second statement, we use again that j = to infer that 

ilK {I-k{v)) ■■= 'tpK = A{iIjk{v)), 

\i£Af / ieAt 

for all V G L^{K\'U3). To prove the third statement, let us consider any g G Pk = 
'ijjx{P)', then using the above definitions we have 

o'if =-i- f {Pi,iJKig))e^Ri)dx = a^{^K{g)) = (TxAg)- 

\K\ Jk 

This proves that P\i{g) = XK{g), hence Ix{g) = □ 

Remark 3.1. The construction of is similar to that used in Girault and Lions 
[isl Appendix] to extend the Scott-Zhang quasi-interpolation operator and make it 
stable in L^{D) for scalar-valued functions. 

3.2. Error estimates for We establish in this section error estimates for 
the operator 

Theorem 3.2 (Local interpolation). There exists a uniform constant c such that 
the following local error estimate holds: 

|u — (3.6) 

for all I G {0:k + 1}, all m G {0:/}, all p G [l,oo], all v G lL^’^(Ar;R"^), and all 
K€Th. 

Proof Let Z G {0:/c -I- 1} and m G {0:/}. 

(1) Let us set Q{w) := w —X^~{w) for all w G The operator Q is well- 

defined since ^ L^{K;W^). Since all the norms are equivalent in P 

and P C ^ W"^’P{K;W), there exists c depending only on Ush and 

K such that ||2i^(w)|lw™.p(lf;R<!) - '^ll^^('^)llLi(i?;R 9 ); which in turns implies that 
ll^^('^)llw™.p(/f-R 9 ) — ''II'^IIli(k-R‘j)’ si'^ce we have already established that is uni¬ 
formly bounded over L^K-W)] hence, G G /:(W'’P(.A; R«); W”"’P(iL; R«)). Assume 
first that I > 1, then [P;_i]'J is point-wise invariant under 1?^ since Z — 1 < fc and 

[Pi-ij'J C [Pfc]'* C P; this in turn implies that the operator G vanishes on [Pi_i]'J. As 
a consequence, we infer that 

|w — 2il'R'lw"*.P(P;R9) = l^('^)lvE™.P(P;R9) “ pg[ipj^i]9 P'^^W^’P{K-,S.9) 

— ll^ll£(W‘.P(if;R'!);VE™'P(P;R'!)) VE'-p(P;R‘i) 

for all w G W^’P{K;'Rg), where the last estimate is a consequence of the Bramble- 
Hilbert/Deny-Lions Lemma. Finally, the above inequality is trivial if Z = m = 0. 

(2) Now let V G W’''P{K-,MP). Using the above argument together with the fact that 
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iff commutes with (see Proposition 13.IL we have 

I'V — ^kV\w^’P{K iRi) < li’K \c{W’P'-P(K;R’i);W’P'P{K;R'i)) IV'A'(f) ~ '4’k{Ik'^)\iy'^,p{K-,Ri) 

— \'’Pk I£(W'"’P(A;R9);W™.!>(A;B9)) IV’A'('*^) “ (i')) | lym.p 

— ^\'4’k \c{W’P’P{K-,Ri)-,W‘^’P{K;Ri)) \'^K:('^)\w''’P{K;Rp) 

^ ^\'^K \c{WPP-P{K-Rp)-W^’P{K-,RP))\'^K:\c(W‘’P{K-,RP)-,W‘-PiK-,Rp))\'‘^\w‘’P{K-RP)- 

The estimate (13.611 follows by using (12.1011 . □ 

4. Averaging operator. In this section, we introduce a bounded linear operator 
jav . P{Th) based on averaging. 

4.1. Connectivity array. Let {‘Pa^a^Ah be a basis of P{Th)\ the functions tfa 

are called global shape functions. We assume that this basis is constructed so that 
for any K € 77i, either int(A) fl supp((/3q) = 0 or there is a unique i £ M such that 
fa\K = (Recall that this is the usual way of constructing finite element bases.) 
We denote by a : T^xW Ah the map such that this map is 

henceforth called connectivity array. Note that a is surjective by definition, but in 
general a is not injective. Denoting by a“^(a) the preimage of a £ Ah, we define the 
connectivity set Ca C Th xA/” for any a £ Ah such that 

Ca ■■= a~^(a) = {{K,i) £ Th^M \ a = a(A,i)}. (4.1) 

Remark 4.1. (Particular case card(Ca) = IJ Assume that card(Ca) = 1, i.e., Ca = 
{(Ao,io)}- Let K 7 ^ Kq, then it is not possible to find an index i € J\f such that 
a = a(A, i) since card(Ca) = 1; this in turn implies that int(A") nsupp((/?a) = 0, hence 
‘Ca\mi(K) ~ b. This means that Lpa is supported on one element only, i.e., ipa is the 
zero extension of Okom- Given the characterization of P{Th) assumed in (12.1711 . this 
means that the y-trace of (pa on the interior faces of K is zero. 

For any a £ A/j, we set Jy = 0 if card(Ca) = 1. If card(Ca) > 2 we define 
Pa G Ph to be the set of the interfaces F such that there are {K,i), {K\i') £ Ca so 
that F = dK r\dK'. We now relate the y-traces to the degrees of freedom by making 
the following assumption: there exists a uniform constant c such that the following 
holds for all v in P'°{Th) and all a & Ah such that card(Ca) > 2: 

|crif.i(u|if) - (JK',i'{v\K')\ < c min(||Aif||^ 2 , ||Aif/||| 2 )|||u];^||ioo(p..Rt), (4.2) 

for all F G P° and all pairs {K,i),{K',i') £ Ca such that F = dK n dKC The 
estimate (SJ) is natural since the degrees of freedom of finite elements providing 
some conformity in JT(curl), or iT(div) are devised to coincide across interfaces 
when the y-jump is zero and when the degrees of freedom in question belong to 
the same connectivity class. In particular, the estimate (14.2p holds true for all the 
finite elements considered in >12.41 Owing to (I2.17p . the assumption (14.21) immediately 
implies that 

\crK,i{'>^\K) - <^K',i'(.'^\K')\ =0, Vu £ P{Th), Va £ Ah, y{K,i), {K',i') £ Ca- 

(4.3) 

4.2. Averaging operator. We define the operator : P^[Th) P{Ph) by 


Jh^'iv) = E 

aeAh 


E ^K,i{viK) 'Ca- 


card(Ca) 


(4.4) 
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For any K gTh, we introduce the notation 

Tk ■= G7h\3i' i)}, (4.5a) 

Dx ■= intiUK'eTK^'}- (4.5b) 


The set Tk is the union of all the cells that share global shape functions with K and 
Dk is the interior of the collection of the points composing the cells in Tk ■ 

Lemma 4.1 (L^’-stability). There exists a uniform constant c such that 


\\Jr(.'^)\\LPiK-,R‘i) < C ||t;||LP(Dj^;R,), (4.6) 

for all p £ [l,oo], all v £ P^{Th), and all K £ Th- 

Proof. We prove the result for p = oo; the other cases are obtained by using local 
inverse inequalities in P'°{Th). Using the triangle inequality and the shape-regularity 
of the mesh sequence {Th)h>o, we infer that 


II ^avr...M| / \\0k,i\\l°°{K-Ri) 

\\Jh (^^) < >,-T 77 ;-^ 

^ card(Ca(if,q) 


E 

{K' 


CTif'.r' {V\K')\ 


^K 


< c ■ 

^^card(Ca(ic.q)^^,^.; 


E \^ K ', i ' iv \ K ')\ 

(A'.»')GC3(K,q 


K'gTk i'GTV 


where the last estimate results from \<Jk' ,i'{'V\K')\ ^ \’^k' ,i'{v\K') ~ (^K,i{v\K)\ + 
WkA'^\k)V the as sumpt ion the inequality |H< IHf|l~(E;R*) < 

lkllL“(DK;R«) and ( |2.12p . □ 

Lemma 4.2 (Approximation by averaging). There exists a uniform constant c 
such that the following holds with the notation := ^y 


1 ?^ — Jh^{v)\wp^’P{K-,Ri) < ch 


If , 1 

a ( — — — + — —m 
Vp r J r 


K 




IIH 


L'-(F;R‘) 


(4.7) 


for all m £ {0: fe + 1}, all p,r € [1, 00 ], all v £ P'°{Th), and all K € Th- 

Proof We only prove the bound for m = 0 and p = r = 00 , the other cases follow 
by invoking standard inverse inequalities. Let v £ P^{Th), set e = v — J^''(v) and 
observe that e £ P'°{Th). Let K €Th. The bound (12.111) implies that 


l|e||L~(A;R9) < c||A ^^||£2 \aK,i{e\K)\- 

ieAf 

Owing to the definition of Jf'' in (14.41) . we first observe that 

°»''-<°l'^>= card(C.«.,) S (4.8) 

4 ’ ^ {K' ,l')^C,(K,i) 

Note that (JK,i{e.\K) = 0 if card(Ca(R:_i)) = 1 (see Remark l4T|) . Let us now consider 
the case card(Ca(R:_i)) > 2. For all {K',i') £ C^(^K,i), there is a path of mesh cells 
in Tk linking K to K' so that any two consecutive mesh cells in the path share 
a common face F £ each face is crossed only once. Furthermore, if 
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{Kr,ir) S Ca(x,i) are such that dKi fl dK-r = F G -y then (j4.2p implies 

that 


WKi,ii{v\Ki) - CrKr,ir{'^\Kr)\ < C min( 111 1 £2 , ||Aif^||f2)|||w]^||j;,oo(p..Rt). 

As a result, 


l|e||L=<=(R-.R,) < c ||AR-^||f2||Ai^||f2 

whence the estimate (I4J1) readily follows since F]^ := UiGA-^a°(A i)’ card(A/’) is uni¬ 
formly bounded, and the mesh sequence {Th)h>o is shape-regular. □ 

5. Quasi-interpolation operator. Let ll : L^{D;R‘^) P^{Th) be such that 
for all AT £ 7^. We now construct the global quasi-interpolation 
operator : L^{D;W^) —>■ P{Th) by setting 

ir:=drroi“. (5.1) 

Note that P{Ph) is point-wise invariant under since P{Fh) is point-wise invariant 
under and Hence, X^ is a projection, i.e., {X^)'^ 

Lemma 5.1 (W™’^-stability). There exists a uniform constant c such that 

\^h^{v)\w^-P{K-S.i) < c\v\wp>.p{Dk-,S.i)^ (5-2) 

for all p £ [1, oo], all m G {0: fc + 1}, all v G and all K G Th- 

Proof. For m = 0, the estimate follows by combining Proposition 13.11 with 
Lemma l4.11 For m > 1, the triangle inequality implies that 

\FT{v)\wp^’P(K-m) < |2ii-('c)|w'".j>(A;R‘i) + \I-\iiv) - Jh"{x\^{v))\w^,p{K-m)- 

Let Ti and T 2 be the two terms on the right-hand side of the above inequality. Ti 
is estimated using Theorem 13.21 with I = m, yielding |Ti| < c\v\wp^.p{K-,s.i)- ^2 is 
estimated using Lemma |42] and the fact that v G hF’"’P(A>R-; R'*) C VF^’^(I1r';R'?) 
has zero 7 -jumps across interfaces (see (12.161) 1. More precisely, we have 

h^\X^\<cl^^ ^ IIPlc(^)l^lUp(F;Rq=c4 ^ II ^ - 4(^)l^lUp(F;Rq 

\ 

^ IK^"-2'^(i^))|A'|Up(F;R‘!) < c/l^ ^ \v\w^.p{K'-,MP), 

K'gTk Fc9A'nF^ F'gTk 

where we have used the triangle inequality to bound the jump by the values over 
the two adjacent mesh cells, the trace inequality from Lemma 17.21 with s = 1, the 
approximation result of Theorem 13.21 and the shape regularity of the mesh sequence. 
Combining the bounds on Ti and T 2 gives (lOI) . □ 

Let us now estimate the approximation properties of Xff'. Since we are going to 
establish error estimates in fractional Sobolev spaces, in the rest of the paper all the 
results that are stated in fractional Sobolev spaces assume that p G [l,oo), whereas 
the results with integer degree of smoothness hold for p G [l,c)o]. Assuming that 
r ^ N and denoting by [rj the largest integer less than or equal to r, we consider 
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the so-called Sobolev-Slobodeckij norm defined as ||r^||vK'’'!>(D;R9) = (II^II^Lrj,p(£).Rq) + 


|r^|vyp,p(_D;R<i) 



Wv{x) - 


||a; 


y\\ 


ir— [rj)p-|-d 
^ 2 (R‘^) 



jl 

P 


(5.3) 


Theorem 5.2 (Local approximation in fractional Sobolev spaces). There exists 
a uniform constant c such that 


\'V "‘|r!|w'-'!>(D*:;R«)! (5-4) 

for all r € [0,fc-|-l], all m G {0: [rj}, all p G [l,oo) ifr^N or all p G [l,oo] i/r G N, 
all V G lT’'’P(i:>R-;R«), and all K G %■ 

Proof Using that Iff {g) = g for all g G together with the stability of Iff 

in the lU'”’^-seminorm and the triangle inequality, we infer that 

- Iff = \v - g-Iff (v - g)\wm,p(K-;Ri) 

< k — ff|vE'".P(if;R‘!) + \Iff{'^ — g)\w^’-p{K--Ri) 

<c\v — 5 |w'".J>(£>k-;R‘>)' 

That is to say, \v-Iff (v)\wm.P(K-m.‘i) < cinfge[Pj, \v - g\w^^p(DK\M.‘i)- We conclude 
by applying Lemma 15.31 componentwise. □ 

Lemma 5.3 (Polynomial approximation in Dk)- The following holds: 

inf \v - glw^p.piDK) <chff^\v\w-.p{DK)^ (5-5) 

9^^k,d 

for allr G [0, /c-l- 1], all m G {0: [rj}, all p G [1, oo) i/r ^ N or all p G [1, oo] if r G N, 
all V G W^’PiDx), and all K gTh- 

Proo f. We proceed as in Bramble and Hilbert [3, Thm. 1], but instead of invoking 
Morrey [l9l Thm. 3.6.11], where the constants may depend on Dk, we are going 
to track the constants to make sure that they are independent of Dx- If to = r, 
there is nothing to prove. Let us assume that m < r. Let £ G N be such that 
£ = r — lifrisa natural number or £ = [rJ otherwise (note that 1 < r if r is a 
natural number since we assumed that 0 < to < r). In both cases the integer £ is 
such that m < i < k. Let Ae^d = {a G j |Qf| := ai + ... + ad < £}. Note that 
card(A,d) = dim(P^,d) = =: Ng^d- Since the mapping ^g^d ■ ^e,d -t such 

that $Gd(g) = iJjj^d^'qdx) aeAe d is isomorphism, there is a unique polynomial 
T^i{v) GFg^d such that ‘^g,d{T<'e{v)) = (/^^ 5“z;da::)aeA«,d> i-e -,/dk ^ ® 

for all a G Ag^d (this result is actually stated in Morrey [3 Thm. 3.6.10]). 

Since by definition d°‘(v — 7rg(v)) dx = 0 for all |a| = to < £, we can apply 
Lemma 15.51 below, i.e., there is a uniform constant c such that ju — 7:g{v)\w^,p{DK) — 
chxlv — Trg(v)\]jym+i,p(^jgi^y Wc can repeat the argument if to-|- 1 < £ since in this case 
we also have 9“(z; — Trg{v)) da: = 0 for all jaj = to -|- 1 < £. Eventually, we obtain 

k - < c/l^™|u - Trg{v)\wi.p(^DK)- 

If r is a natural number, then £ -|- 1 = r, and we can apply the above argument 
one last time since d°‘{v — Trg{v))dx = 0 for all ja] = £, which then gives (15.5p 
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because d°‘TTf^{v) = 0 for all \a\ = 1+ 1. Otherwise, ^ = [rj and we apply Lemma [7?T] 

to all the partial derivatives d°‘{v — 7r^(u)) with \a\ = £, s = r — [rJ S (0,1) and 

O := Dk', this is legitimate since all these partial derivatives have zero average over 

Dr- We infer that there is c uniform with respect to s, p, K, and v such that 

1 
P 

|u — 7r^(u)|w'-.p(nK;R'')- 

Note that \v— Tri{v)\w’'.p{DK■,«■’’) — klw'’'P(DK;R‘*) since d°‘'!re{v) is a constant in R'J for 
all |a| = £. We conclude that (j5.5ll holds owing to the shape-regularity of the mesh 
sequence. □ 

Corollary 5.4 (Global best approximation in L^). There exists a uniform 
constant c such that 


|u — 7rf(u)|vym,p(£,^.R,) < ch 


[rj-m, r-[rj 


K 


''Dk 


hloK 

\Dk\ 


inf \\v - WhWLPiD-^Ri) < ch''\v\wp.p(D-,Ri), (5.6) 

wh&PiTh) 

for all r G [0,fc -|- 1], all p € [1,cxd) if r ^ N or all p G [1,oo] if r G N, and all 
V G W''’P{D-,W). 

Remark 5.1. (Interpolation) Corollarv l5.4l can also be proved using (15.41) for r G N 
and the real interpolation method (i.e., the JlT-method), see e.g., Tartar [2J, Chap. 22]. 

Remark 5.2. (Approximation for Note in passing that Theorem 13.21 which 
states the approximation properties of , can be re-written with fractional Sobolev 
norms, i.e., the following also holds: 

\v — ^r{v)\wpP’P{K;R'i) < (5.7) 

for all r G [0, fc -I- 1], all m G {0: [rJ}, all p G [1, oo) if r ^ N or all p G [1, oo] if r G N, 
all V G W^’P{K;W), and all K gTh. 

Lemma 5.5 (Poincare inequality in £)if). Letwjj^ be the average ofw over Dr. 
There exists a uniform constant c such that 

h-V r^Wlp^Dk) <chR\v\wrp(DK)^ ( 5 - 8 ) 

for all p G [1, oo], all v G W^’P(Dr), and all K gTh. 

Proof Let K G Th- Let Ki,Kr G Tr sharing an interface F = dKi n dKp. We 
observe that 

\Vr, -VrJ = \F\-p\\Vr^ -U|if, +Vir^-VrJ\lp(F), 

since v^r^ = v\r^ on F. By using the triangle inequality, estimating the two norms 
in LP{F) with the trace inequality (17.211 (with s = 1), and by applying the Poincare 

inequality in Ki and Kr separately (both with constant Tr~^ since the mesh cells are 

1_1 1 — — 

convex sets), we obtain \vr^ — VrJ < c(/i^^ ” |i;|wi,p(if,) + Iir^^\ v\w^’P{Kr))- After 
invoking the shape regularity of the mesh sequence, we infer that 

\K\p\vr^ -VrJ < chR\v\wrp(RiuKr.)- (5-9) 

Observing that - Vr, = Y.r„^Tk ~ l^K') for any K' G Tr, we infer 

that 

h-Vjj^\\LP{R')<\\v-VR,\\LP(R')+ ^ T7^,\K'\p\Vr„ -Vr,\. 

K"£Tk ' 
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For any K" G Tk, we can find a path of mesh cells in Tk linking K' to K" so that 
any consecutive mesh cells in the path share a common face and this face is crossed 
only once. Using (ICTD together with the shape regularity of the mesh sequence, we 
infer that ||u — ‘}Ik’\\l'p{k') ^ and the conclusion follows by summing 

over K' gTh and using the fact that card(TR-) is uniformly bounded. □ 

6. Quasi-interpolation with bonndary prescription. Our goal in this sec¬ 
tion is to construct a variant of the quasi-interpolation operator that prescribes 
homogeneous boundary values. 

6.1. Trace operator. Let F G be a boundary face. We denote hy Kp 
the unique cell such that F C dKp. We consider the global trace operator 7 : 


-G L^{dD-,R*) such that 

liv)\F =lKp,Fiv\KF), VUeJ'f. (6.1) 

We define for all p G [1, c)o] the functional space 

:= {v G \ y{v) = 0}. (6.2) 

We then set 

PoiTh) ■■= {vh G P{Th) I l{vh) = 0}. (6.3) 

The typical examples we have in mind are 

Pl{Th) := {vh € P^{Th) I = 0 }, (6.4a) 

P^iTh) := {vh G P^{Th) I = 0 }, (6.4b) 

PoiPh) := {vh G P'^iTh) I Vh-n\gp, = 0}. (6.4c) 

Upon setting V)f = {w G U® | v\gp, = 0}, = {u G | vxn^gp, = 0}, and 


g yd I ^ 0 }, we have P^iTh) = PHTh)nVo^, P§ijl) = P^{Th)r^VS, 

and PoiTh) = P’^iTh) r\ i.e., the finite element spaces {%), P§ {%), Pq {%) 
are conforming in Vq, , respectively. 

In the rest of the paper, we slightly abuse the terminology by calling global degrees 
of freedom the elements of Ah- We say that a global degree of freedom a G Ah '^s an 
internal degree of freedom if yi^Pa) = 0. The collection of all the internal degrees of 
freedom is denoted .A)); the degrees of freedom in Af = Ah\Af^ are called boundary 
degrees of freedom. Let a G Af^, then there is a face F G such that 'y{ipa)\F T 0- 
Let Kf be the unique cell such that F C 9/Ff, then supp((/Jq) r\KF ^ 0. This means 
that there is a unique if G Ff such that iPa\Kp = (^Kp,i- For all a G A^, we define 
to be the collection of all the boundary faces F G such that there is [KfAf) G Ca 
and F C dKp] we set J'f = 0 if a e A.)). We abuse the notation by setting 

Mpix) =yKp,Fiv\Kp)ix), and ^^(a;) = U|F(ic), as a; € F, VU G J'f, (6.5) 

and assume that ||u]^(a;)| < ||w]F(a;)|, sex G F, for all F G T^. In coherence with the 
assumption (03), we assume that there is a uniform constant c such that the following 
holds for all the boundary degrees of freedom a G Af, all F G Ff, alH^ G W such 
that (KpAf) G Ca, and all v G P'°{Th)- 

WKp,ipiv)\ < c ||AFF||f2||7F'p_F(f|Fp.)|U“(F;R‘), (6.6) 

Note that this assumption is satisfied by all the finite elements considered in H2.4I 
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6.2. Averaging and quasi-interpolation operators revisited. We are going 
to modify the averaging operator to prescribe homogeneous boundary conditions. 
We define ■ P'°{Th) ^ Po(Jh) by setting for all v G P'^{Th), 


Jhoiv) = E 


aGA' 


; 1 card(C, 


E ^kA'^Ik) ‘Pa- 
{K,l)GCa. / 


Lemma 6.1 (L^-stability). There exists a uniform constant c such that 


(6.7) 


\\>^hoi'‘^)\\LPiK-,Ri) < C ||u||iP(£)^.R 5 ), 


( 6 . 8 ) 


for all p G [l,C)o], all v G P'°{Th), and all K G Th- 
Proof. Proceed as in the proof of Lemma 14.11 □ 

Lemma 6.2 (Approximation by averaging). There exists a uniform constant c 
such that the following holds: 


\w 


< ch 


\p 


K 




(6.9) 


for all m G {0:k + 1}, all p,r G [l,oo], all w G P'°{Ph), and all K G Th, with 
Tk — and defined in Lemma 14-^ 

Proof. This is a straightforward adaptation of the proof of Lemma 14.21 Letting 
e = V — (u), we observe that (jK,i{e\K) is still given by (14.81) if a(A, i) G A°f^, while 

<yK,i{a\K) = o'K.iiw) if a{K,i) G Af, and this term is bounded using dSH). □ 

A global quasi-interpolation operator : L^(H;R'^) —>■ Po{Th) is then dehned 
by setting 

= ( 6 - 10 ) 

Note that Po{Th) is point-wise invariant under I))q since ()6.9() implies that Po{Th) is 
point-wise invariant under Hence, is a projection, i.e., =^ho- 

Lemma 6.3 (L^’-stability of I^q). There is a uniform constant c such that 

\\'^hoi'a)\\LP{K-Ri) < C ||u||lp(Dj^;R<!), (6.11) 

for all p G [1, oo], all v G L^{D] M"^), and all K G Th. 

Proof. Proceed as for the proof of Lemma 15.11 □ 

6.3. Error estimates. The purpose of this section is to establish error estimates 
for the quasi-interpolation operator I))g in the W’^-norm (either integer or fractional). 
Let r G [0, fc-l-1] and p G [1, oo) if r ^ N or p G [1, oo] if r G N. If r > i, then functions 
in IV’''^(T);R'^) have traces on dL>, and therefore it makes sense to define 


:= {v G W^’^iD-R^^) \ -fiv) = 0}. (6.12) 

We are going to use the following notation 

T,° := {A G I Vi G AT, a(A,i) G AU, (6.13a) 

Th^ ■■= ThW = {KGTh\3iGN, a{K,i) G At}, (6.13b) 

:= int(U;fg 7 -sA). (6.13c) 
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7^° is the set of the cells whose degrees of freedom are all internal. is the set of 
the cells that have at least one boundary degree of freedom. is the collection of 
the points in D that belong to at least one cell in . 

Theorem 6.4 (Approximation). There exists a uniform constant c such that the 
following estimate holds for all r £ [0, fc + 1], all m G {0: [rj}, all p £ [1, oo) if r ^ N 
or all p £ [l,oo] j/r £ N, all v £ and all K 

1^ “ (6-14) 

Moreover, (j6.14p also holds if rp > 1 with c depending on \rp — 1| for all v £ 
Wq’^{D;W^) and all K £ 7)f. Finally, if rp < 1 (i.e., r £ (0,1) and m = 0), we 
have 

\\v-^hoi'^)\\LP{Da-,S.<i) < ch^\\v\\w-.P^D;Ri), G W^’P{D;W). (6.15) 


Proof Let AT be a cell in Th- If AT £ then I^q{v)\x = this 

proves (16.1411 in this case. Let us now consider K G 7)f. The triangle inequality 
implies that 

l^hoM - < l^h''('a) - v\w^.p{K-,m) + \^ho{v) - i‘v)\wp‘^p(K■,R1)■ 

Sm.ce we have already established that \I^(v) — v\wm,p{K;R’i) < ch^x^"^\v\wp,p{DK-,s.i) 
in Theorem O we just need to estimate — Iff {v)\wp^<p{k-;k.'i)- Let us define 

the set of the boundary degrees of freedom with nonempty support on AT, := {a £ 
Af I £ A/", a(Ar, i) G Ca}- Then 


(ifiiv) -xr(d)iic = - 






Ok,, 


For any a G A% and any (AT', i') G Ca, there is a face F G A’f and a pair {Kp,iF) G Ca 
such that there is a path of mesh cells in Tk linking K' to Kp so that any two 
consecutive mesh cells in the path share a common face in iF°, and each face is 
crossed only once. This observation implies that 


k(if/,*/)(g)| < 


Wa(Ki,ii)i9\Ki) - CPB{Kr,ir)i 9 \Kj\ + ka(ifF .ip) (ffliCp ) I 


for all g G P^{7h)- By proceeding as in the proof of Lemma [4.21 we infer that 


+ X! II'7(2^*(^^))|f||lp(F;R*) ■ 

F&FS 

Case 1, rp > 1: Let us assume that v G The boundary condition 

'y{v)p = 0, for all F G and the continuity condition |n]J. = 0, for all F G Fff, 
imply that 

“2:,7(u)|wm,p(^.R,) < ^ II(•(;)-■!;||iP(F.Rt). 

FeF^UF® 
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The conclusion follows by invoking the trace inequality from Lemma 17.21 either with 
s = 1 if r € N or with s = r — [rj otherwise, and the approximation properties of 
stated either in Theorem 13.21 or in Remark |5.21 

Case 2, rp < 1: Assume now that rp < 1. Norm equivalence implies that 

Let p be the distance to dD] then there is c uniform with respect to the mesh sequence 
such that ||p||j;,oo(£ia) < ch and 

l|2iho(^) ~ ^h'i'^)\\LP{D^-,m) < c(||X^(u) — r'||LP(Da;R<z) + ||l'||LP(Da;R9)) 

< C + ||p’'p ’'^^||LP(Da;R<i)) 

< c ^/l''||u||wp.p(D;R‘!) + l|p|lL=x=(£) 8 )||p“"'?^||LP(Da;R 5 )^ . 
Since rp < 1, we infer that (see Grisvard [3 Cor. 1.4.4.5]) 

Up ’''<^llLP(Da;R9) < Up '’^^||lp(D;R 9) < c||u||vyp,p(£).R9). 

In conclusion, ||X^^(r;) - Tri'^)\\LPiDS-mi) < ch^\\v\\w'--piD-,m)- □ 

Corollary 6.5 (Global best approximation in L^). There exists a uniform 
constant c, additionally depending on \rp — 1|, such that 


inf ||u - u;/i||z,p(r(;R‘!) < 
’Wh&PoiTh) 


c/l’’|u|n/P,P(_D;R9), 

ch^\\v\\wp.P{D-R'i), 


V-y G Wq’^{D;W) ifrp > 1 
Vw G 1T’'’P(X»;IR«) ifrp < L 


Remark 6.1. (Theorem \6.4\^ For rp > 1, a similar estimate has been obtained in 
the scalar-valued case for the Scott-Zhang interpolation operator by Ciarlet Q . Fur¬ 
thermore, the estimate for rp < 1 in Theorem 16.41 essentially says that the difference 
V — T^q{v) does not blow up too fast close to the boundary. A better result is not 
expected since T^q{v) is forced to be zero at dD whereas v can blow up like p~‘‘w 
where w is a function in LP{D;R.‘^). 

Remark 6.2. (rp = 1) Let r G (0,1). Using the notation from the real interpola¬ 
tion th eory , it is known that W^’P{D) = [X^(X>), W^’P{D)]r^p since D is Lipschitz, see 
Tartar 24, Lem. 36.1]. Let us define 

lUo’’(;(i^;M^) := [XJ’(U; K«), (6.17) 


Then, using Theorem 16.41 with I G {0,1} and m = 0, the real interpolation theory 
implies that 


11 ^ “ 2 i))o('y)||iP(£|.R<,) < C/l’'||y||iyPj;P^(D;Rq), (6.18) 

for all p G [1, cx)) and all v G R**). This estimate is not fully satisfactory for 

two reasons. First it is not local. Second it is not really clear what Wqq ,^,(1?; R'^) is. 
For instance, let us define 

lFo'’^’(i^) := {v G W^’P{D) I v^ 9 D = 0}, (6.19a) 

W^’^{D) := {u G W^'P{D) I vxTUgo = 0}, (6.19b) 

W}f^{D) := (u G W^^P{D) I vruoD = 0}. (6.19c) 

One then realizes that characterizing [LP{D),W^’^{D)]i p and [LP{D),W^^{D)]i p 
in terms of Sobolev regularity is (possible but) not straightforward, and to the best 
of our knowledge, a full characterization of these spaces is not yet available. 
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7. Technical results in fractional Sobolev spaces. This section contains 
two technical results in fractional Sobolev spaces: a Poincare inequality and a trace 
inequality. 

Lemma 7.1 (Poincare inequality). Let O be an open set in and let Vq be 
the average of v over O, for any v £ L^(0). Let ho '■= diani(O). Then, for all 
V £ W^’P(O) with s £ (0,1) and p £ [1, oo), the following holds: 

/ \ p 

Ik ~ LoIIl!’( 0) < klvE'>’P(0)- (7.1) 


Proof. This result is essentially Proposition 6.1 from Dupont and Scott [l^, see 
also Heuer 17|. We nevertheless give a proof since the computation gives us the 
constant in the right-hand side of (EU, and this in turn allows us to apply the result 
locally on shape-regular meshes. Using the definitions, we have 

f \v{x)-Vofdx= [ \0\~P [ {v{x)-v{y))dy 

Jo Jo Jo 

k(ic) - v{y)\ 


dx 


< / ior" 

Jo 


Jo 




x-y\ 


nS+-, 


< [ |o| 

Jo 


-p 


x-y\ 
r \v{x) - v{y)\P 
lo \\x-yfP'^'‘ 


dy 


dy 

j da: 


is+§)p' 

y|| 



dy dx, 


where p' := Then using that ||a; — y \\^2 < ho for all x,y € O, we infer that 


Ik-Lollip(o) ^ / 1^1 ^ 


r \v{x) - v{y)\v 
'o \\x-yf^^+‘^ 


dydximax/ ||a; - y ||^2 


(s+|)p' 


yajGO Jq 


dy 


<krw»,.(0)|OrM I ho ^ dy 


, (s-hI)p' 




Hence Ik-^olUuo) < (|^) " \v\w‘’P{ 0 )- D 

Lemma 7.2 (Trace inequality). Assume s £ (0,1) and sp > 1 with p £ [l,oo) or 
s = 1 with p £ [l,oo]. There exists c, uniform with respect to the mesh sequence but 
depending on \sp — 1| if s G (0,1), such that the following holds for all v £ 
and all K gTh: 

IkllLp(E) < c(/l^*’ |k||_LP(if) + hj^ ” \v\w=,p(K))- (7.2) 


Proof Let v £ W^’^{D). Let K € Th he a cell and F be a face of K. Since 
the map Tk is affine, using a trace inequality in W^’P{K) (recall that s £ (0,1) and 
sp > 1 or s = 1 and p > 1), we infer that 

I F I p 1 

IklUu^’) = 7^llki(v)lkp(P) < Cs.p|F|p(||V’|(z;)lkp(jf) -f |V'|(^^)lw'=.P(A))’ 

|F|p 
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where Cg^p depends on |sp — 1| if s S (0,1). Upon changing variables, this inequality 
is re-written 


\F\-.\K\-H\\v\\ LP{K) + pA'||f2*k|rV'»>P(if))- 

The conclusion follows from the shape-regularity of the mesh sequence, see (1^ . □ 
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